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Abstract — We investigate the problem of a monostatic radar 
transceiver trying to detect targets, sparsely populated in the 
radar's unambiguous time-frequency region. Several past works 
employ compressed sensing (CS) algorithms to this type of 
problem, but either do not address sample rate reduction, impose 
constraints on the radar transmitter, or propose CS recovery 
methods with prohibitive dictionary size. Here, using the Xam- 
pling framework, we describe a sub-Nyquist sampling approach 
which has several important advantages over previous methods: 
low rate sampling and digitai processing, undersampling in both 
fast and slow time, the resulting CS dictionary size is proportional 
to delay grid size only, and we impose no restrictions on the 
transmitter. Xampling allows reducing the number of samples 
needed to accurately represent the signal, directly in the analog- 
to-digital conversion process. Furthermore, our approach can be 
implemented in hardware using a previously suggested Xampling 
prototype. After sampling, the entire digitai recovery process is 
performed on the low rate samples without having to return 
to the Nyquist rate. In the noiseless case, we can recover a 
sparse target scene from samples taken at a rate proportional 
to the signal's rate of innovation, which can be several orders 
of magnitude smaller than its Nyquist rate. In the presence of 
noise, we introduce the concept of Doppler focusing from the low 
rate samples in order to obtain good detection performance even 
at SNRs as low as -25dB. 

Index Terms — compressed sensing, rate of innovation, radar, 
sparse recovery, sub-Nyquist sampling, delay-Doppler estimation. 



I. INTRODUCTION 

We consider target detection and feature extraction in a 
radar system, using sub-Nyquist sampling rates. The radar is a 
single transceiver, monostatic, narrow-band pulse-train system. 
Targets are non-fluctuating point targets, sparsely populated in 
the radar's unambiguous time-frequency region: delays up to 
the Pulse Repetition Interval (PRI), and Doppler frequencies 
up to its reciprocai the Pulse Repetition Frequency (PRF). 
We propose a recovery method which can detect and estimate 
targets' time delay and Doppler frequency, using a linear, non- 
adaptive sampling technique at a rate significantly lower than 
the radar signal's Nyquist frequency, assuming the number of 
targets L is small. 

Current state-of-the-art radar systems sample at the signal's 
Nyquist rate, which can be hundreds of MHz and even up to 
several GHz. Systems exploiting sub-Nyquist sampling rates 
can benefit from a lower rate analog-to-digital conversion 
(ADC), which requires less computational power. Moreover, 
sampling at the Nyquist rate may not always be feasible due 
to high power consumption, heat dissipation, cost, or other 
practical considerations. Finally, offiine radar systems which 
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record samples for subsequent processing will gain substantial 
Storage capacity reduction if they are able to sample at sub- 
Nyquist rates. 

The goal of this work is to break the link between radar 
signal bandwidth and sampling rate. The sub-Nyquist Xam- 
pling ("compressed sampling") JT] method we use is an ADC 
which performs analog prefìltering of the signal before taking 
point-wise samples. These compressed samples ("Xamples") 
contain the information needed to recover the desired signal 
parameters using compressed sensing (CS) algorithms. This 
work expands the work in 0, adding Doppler to the target 
model and proposing a new method to estimate it. The same 
sampling technique and hardware that were used in Q, O 
will also work here, while the digitai processing we suggest 
is different. 

Past works employ compressed sensing (CS) algorithms to 
this type of problem, but do not address sample rate reduction 
and continue sampling at the Nyquist rate H, 0. Other 
works combine radar and CS in order to reduce the receiver's 
sampling rate, but in doing so impose constraints on the radar 
transmitter and do not treat noise (6), or require an infinite 
number of samples [7]. Another line of work proposes single 
stage CS recovery methods with dictionary size proportional 
to the product of delay and Doppler grid sizes, making them 
infeasible for many realistic scenarios |5|, |8|. 

Our approach is based on the observation that the received 
radar signal can be modeled with 3L degrees of freedom 
(DOF): a delay, Doppler frequency and amplitude for each 
of the L targets. Signals which can be described with a fixed 
number of DOF per unit of time are known as Finite Rate of 
Innovation (FRI) IH signals. The proposed recovery process 
is actually a recovery of these DOF from low rate samples. 
The concept of FRI together with the Xampling methodology 
enables sub-Nyquist rates using practical hardware (TJ. 

At the crux of our proposed method is a coherent super- 
position of time shifted and modulated pulses, the Doppler 
focusing function For any Doppler frequency v, this 

function combines the received signals from different pulses 
so targets with appropriate Doppler frequencies come together 
in phase. For each sought after v, §(t;v) is processed as a 
simple one-dimensional CS problem and the appropriate time 
delays are recovered. The gain from this method is both in 
terms of signal-to-noise ratio (SNR) and Doppler resolution. 
For P pulses adding coherently, we obtain a factor P SNR 
improvement over white noise (which adds incoherently, i.e. 
in power). In addition, denoting the PRI as r, the width of the 
Doppler focus for each $(£; v) is 2-k/Pt, meaning that delays 
of targets separated in Doppler by more than 2tt/Pt will not 
interfere with each other. 

The idea for Doppler focusing comes from a similar func- 
tion in the context of ultrasound beamforming used in |10|. 
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There, in a method named "Dynamic focusing", the signal 
returned to a set of linearly aligned transceivers is focused in 
a manner similar to how we focus pulses, and the Doppler 
frequency v is replaced by spatial direction 9. In both cases, 
the advantages of focusing are not lost with sub-Nyquist 
processing since it can be performed on the low rate Xamples 
and improves the SNR. 

Simulations provided in Section |VII| show that when sam- 
pling at one tenth the Nyquist rate, our Doppler focusing 
recovery method outperforms both two-stage CS recovery 
and classic processing. When the SNR reaches -25dB, our 
approach achieves the performance of classic matched filtering 
(MF) operating at the full Nyquist rate. 

The main merits of our proposed method are as follows: 
1) Low rate ADC and DSP - using Xampling and the 
proposed recovery method, we are able to acquire the 
sub-Nyquist samples containing information needed for 
target recovery, and then digitally recover the unknown 
target parameters using low rate processing, without 
returning to the higher Nyquist rate. 
Undersampling in both fast and slow time - our 
method allows for sample rate reduction both in the 
intra-pulse ("fast time") and the inter-pulse ("slow 
time") samples. This provides an additional degree of 
freedom in the algorithm, facilitating improved design of 
system architecture and hardware. For example, systems 
with a constraint on average (as opposed to instanta- 
neous) sample rate may choose to alleviate the fast time 
sample rate reduction at the expense of reducing slow 
time sample rate. 

Scaling with problem size - many CS delay-Doppler 
estimation methods depend upon constructing a CS 
dictionary with a row for each delay-Doppler hypothesis. 
For even moderate size problems, the number of delay 
or Doppler hypotheses (i.e. grid sizes) can easily be 
on the order of IO 3 , with a number of measurements 
of similar order. This entails the dictionary to posses 
IO 9 elements, which is infeasible for many systems, 
especially real-time ones. Our Doppler focusing based 
method avoids this problem by separating the Doppler 
from delay recovery, making each CS delay recovery 
indifferent to the underlying Doppler. 
Transmitter compatibility - our recovery method does 
not impose any restrictions on the transmitted signal, 
provided it meets the assumptions stated in Section [II] 
The remainder of this paper is organized as follows. In 
Section [II] we describe the radar model and the assumptions 
used for simplifìcation. Section [III] reviews classic matched 
filter processing. We explain the Doppler focusing concept 
in Section IV and sub-Nyquist delay recovery in Section [V] 
The delay-Doppler recovery method using Doppler focusing 
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is described in Section VI and numerical results are presented 
in Section IVIII 

We denote vectors by boldface lower case letters, e.g. c, 
and matrices by boldface capital letters, e.g. A. We say a 
vector c is L-sparse if ||c|| < L, i.e. at most L of its 
elements are nonzero. The nth element of a vector is written 
as c„, and the ijth element of a matrix is denoted by Ay. 



Non-boldface variables represent scalars or functions, where 
continuous functions are denoted with round parentheses, e.g. 
x(t) and discrete functions with square parentheses, e.g. h[n]. 

II. Radar Model 
We consider a radar transceiver that transmits a pulse train 

p-i 

x T (t) = ^2h(t-pT), 0<t<Pr (1) 

consisting of P equally spaced pulses h(t). The pulse-to- 
pulse delay r is referred to as the PRI, and its reciprocai 
1/r is the PRF The pulse h(t) is a known time-limited 
baseband function with continuous-time Fourier transform 
(CTFT) H(uj) = h{t)e- iut db. We assume that H(w) has 
negligible energy at frequencies beyond Bh/2 and we refer to 
P>h as the bandwidth of h(t). The target scene is composed of 
L non-fluctuating point targets (Swerling model, see ifTTI '). 
where we assume that L is known, although this assumption 
can easily be relaxed. The pulses reflect off the L targets and 
propagate back to the transceiver. Each target l is defìned 
by three parameters: a time delay ti, proportional to the 
target's distance from the radar; a Doppler radiai frequency v\, 
proportional to the target-radar closing velocity; and a complex 
amplitude ai e C, proportional to the target's radar cross 
section (RCS), dispersion attenuation and ali other propagation 
factors. We limit ourselves to defìning targets in the radar's 
radiai coordinate system; future work will deal with targets' 
angular properties. 

Throughout, we make the following assumptions on the 
targets' location and motion, which allow us to obtain a 
simplifìed expression for the received signal. See the Appendix 
for details. 

Al. "Far targets" - target-radar distance is large compared 

to the distance change during observation interval which 

allows for Constant ai. 
A2. "Slow targets" - small target velocity allows for Constant 

ti during observation interval and Constant Doppler 

phase during pulse time T p . 
A3. "Small acceleration" - target velocity remains approxi- 

mately Constant during observation interval allowing for 

Constant v\. 

Although these assumptions may seem hard to comply with, 
they ali rely on slow "enough" relative motion between the 
radar and its targets. Radar systems cracking people, ground 
vehicles and sea vessels usually comply quite easily. For 
example, consider a P=100 pulse radar system with PRI 
r=100^sec, pulse width T p =l/isec, bandwidth _B/ l =30MHz 
and carrier frequency / c =3GHz, tracking cars traveling up to 
120km/hour. The maximal distance change over observation 
interval is approximately 0.33m, so if the targets' minimal 
distance from the radar is a few meters, then Al. is satisfìed. 
As for A2., the maximal Doppler frequency is approximately 
667Hz, which is much smaller than both f c / PtBi^IQKHz 
and l/T p =lMHz. An extreme acceleration of 10m/sec 2 would 
cause a velocity change of O.lm/sec over observation interval, 
easily satisfying A3.. As for airborne targets, care must be 
taken to ensure compliance. 
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Using these three assumptions, we can write the received 
signal as 



P-l L-l 



x(t) = ^ - Tl - Vr)e-^ r . (2) 

p=0 1=0 

It will be convenient to express the signal as a sum of single 
frames 

p-i 

x(t) = x p(t) 

p=0 

where 



(3) 



x p (t) 



L-l 

E 

1=0 



aih(t-n - P T)e- ju ' pT . 



(4) 



In reality x(t) will be contaminateci by additive noise. We will 
take this into account in our simulations. 

Since the transmitted signal ([TJ is a finite periodic pulse 
train, it is invariant to the transformation t; — > ti +k\T, vi — > 
vi + 27rfc 2 /r where ki,k 2 € Z, except on its boundaries. 
Specifically, assuming that k\ <C P, 



P-1L-1 



J9=0 1=0 



-j(vi+2Trk 2 /r)pT 



P-l+fci L-l 



E 



aih(t- ti -pr)é 



«7 



p=kx 1=0 
P-l+fei L-l 

E E 

p=fci 1=0 
P-l+fci L-l 

^ ^aj^t-Ti-jn-Je-^, 



(5) 



with = aie^ lklT . Therefore the radar's unambiguous 
time-frequency region, where it can resolve targets with no 
ambiguity, is [0, r] x [0, 2tt/t] respectively. We make the 
following further assumptions on targets' delay and Doppler: 
A4. No time ambiguity: {tv e / C [0, t)}^, 1 where / is a 
continuous time interval in [0, r), so that x p (t) — 0, Vt $5 
[pr, (p+ 1)t]. 
A5. No Doppler ambiguity: {vi e [0, 27r/r)}^7 1 . 
A6. The pairs in the set \ti,v{\\~^ are unique. 

Our goal in this work is to accurately detect the L targets, 
i.e. to estimate the 3L DOF {ai, Ti, v(\f^ Q in using the 
least possible number of digitai samples. 

III. Matched Filter Processing 

Classic radar processing samples and processes the received 
signal at its Nyquist rate Bh using a MF ifTTIl . In modem 
systems the MF operation is performed digitally, and therefore 
requires an ADC capable of sampling at B^, which can be 
hundreds of MHz and even up to several GHz. In order to 
evaluate our sampling and reconstruction method, we compare 
it to the classic radar processing, which in general consists of 
the following stages: 

1) ADC - sample each incoming frame x p (t) at its Nyquist 
rate Bh, equal to h(t)'s bandwidth, creating a; p [n],0 < 
n < N, where N = tB h . 
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2) Matched Alter - for each x p [n], create y p [n) = x p [n] * 
h* [—n], where h[n] is a sampled version of the transmit- 
ted pulse h(t). The time resolution attained in this step is 
1/Bh, corresponding to the width of the autocorrelation 
of the pulse hit). 

Doppler processing - for each discrete time index 
n, perforai a P-point DFT along the "slow time": 

z n [k] - DFT P {y p [n}} = £p=o y P [n] e -^ k / p for 
< k < P. The frequency resolution attained in this 
step is 2tt/Pt, proportional to the inverse of the total 
coherent processing interval. 

Delay-Doppler map - stacking the vectors x n , and 
taking absolute value, we obtain a delay-Doppler map 
Z = abs[z ... Zat-x] e R PxN . 

Peak detection - a heuristic detection process, where 
knowledge of number of targets, target power, clutter 
location, etc. may help discover target positions. For 
example, if we know there are L targets, then we can 
choose the L strongest points in the map. 
Super-resolution - locai interpolation around detected 
peaks helps achieve sub-pixel resolution. 

Classic processing requires sampling the received signal at 
its Nyquist rate Bh, which is inversely proportional to the 
system's time resolution. The required computational power is 
P convolutions of a signal of length N = tB^ and N FFTs of 
length P - both also proportional to Bh- The growing demands 
for improved estimation accuracy and target separation dictate 
an ever growing increase in signal's bandwidth. The goal of 
this work is to break the link between radar signal bandwidth 
and sampling rate, and to allow low rate sampling and 
processing of radar signals, regardless of their bandwidth. 

We achieve this goal by utilizing the combination of ideas 
of FRI and Xampling. Previous works have already used 
these complementary concepts together. The work in [3|,|9] 
creates a mathematical framework for sub-Nyquist sampling 
of pulse streams and defines lower bounds on the sampling rate 
needed for perfect reconstruction. Practical sampling methods 
achieving these bounds are explained in 0, |[T2l . iflOÌ in 
the context of ultrasound and radar, both without Doppler. 
Another work [7| investigates the delay-Doppler estimation 
problem, but assumes an infinite number of received samples 
and recovers the delays and Doppler frequencies in a two- 



stage process. We will show in Section VII that our Doppler 



focusing recovery method is superior to the two-stage method, 
where mistakes in the first stage propagate to the second stage. 
Combining Xampling and Doppler focusing, our 3L DOF 
input signal |2]) will also enjoy the benefits of sub-Nyquist rate 
Xampling and accurate digitai recovery of the target scene, 
effectively breaking the link between signal bandwidth and 
sampling rate. 

IV. Doppler Focusing 

We now introduce and explain the main idea in this paper, 
called Doppler Focusing. This processing technique uses target 
echoes from different pulses to create a single superimposed 
pulse, improving SNR for robustness against noise and im- 
plicitly estimating targets' Doppler frequency in the process. 
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We point out that stages 2) and 3) in classic processing 
can be viewed as delay processing and Doppler processing 
accordingly. Since they are both linear and time-invariant, they 
can be interchanged, performing the DFT before MF. In classic 
processing, performing MF before DFT decreases computation 
latency, so most practical systems carry them out in the noted 
order. However, when using Doppler focusing, the Doppler 
focusing stage must come before delay estimation and the 
order can no longer be reversed. 

When interchanging steps 2) and 3), the DFT is simply a 
discrete equivalent of the following time shift and modulation 
operation on the received signal: 



p-i 

$(t;i/) = ^ Xv {t + pT y^ 

p=0 

P-l L-l 
p=0 1=0 

L-l P-l 

= ^ a t h(t - n) e i{ - v ~ Vl)vr 

1=0 p=0 



(6) 



where we used Q. 

We now analyze the sum of exponents in Q. For any given 
v, targets with Doppler frequency v\ in a band of width 2it/Pt 
around v, i.e. in v)'s "focus zone", will achieve coherent 
integration and an SNR boost of approximately 



p-l 

p=0 



|i/-i/ ( |<2tt/Pt 



P 



(7) 



compared with a single pulse. On the other hand, since the sum 
of P equally spaced points covering the unit circle is generally 
close to zero, targets with v\ not "in focus" will approximately 
cancel out. Thus g{y\v\) = for \v— v\\ > 2tt/Pt. See Fig.fT 
for an example of g(v\vi). Therefore we can approximate |óf 



200^ 
180 
160 - 
140 




Normalized Doppler frequency v 



Fig. 1: Example of g(v\vi) for P — 200 pulses and vi = 0. 
Red line marks "focus zone", i.e. \v\ < 2-k/Pt. Frequencies 
outside focus zone are severely attenuated. 



(8) 



by 

<$>{t-v) = P J2 aMt-rt). 

l:\v-vi\K2tv I 'Pt 

Instead of trying to estimate delay and Doppler together, 
we have reduced our problem to delay only estimation for a 
small range of Doppler frequencies, with increased amplitude 
for improved performance against noise. To emphasize this, 
consider the case of trying to detect and estimate parameters 
for two targets with very closely spaced delays but with differ- 
ent Doppler frequencies (see for example the two helicopters 
in Fig. |2j. Algorithms whose time resolution is coarser than 
the targets' delay separation are likely to encounter various 
problems recovering this target scene, the most likely of 
which is identification of a single target instead of two. With 
Doppler focusing we achieve an extra dimension of potential 
separation, regardless of the underlying recovery algorithm, 
enabling improved recovery performance. Fig. [2] illustrates 
this concept by showing various targets spanning some delay- 
Doppler region. When focusing for some v only targets in 
i/'s focus zone (white region) come into view, while ali other 



targets (red region) disappear. In Section VII we demonstrate 
this point via simulation. 
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Fig. 2: Schematic delay-Doppler map. Red area indicates out- 
of-focus region. Only one target is in focus for current choice 
of v. 

To allow our sub-Nyquist recovery method, in the next 
sections we sample our signal in the time domain but extract 
frequency domain information. We now show how Doppler 
focusing can also be performed in the frequency domain, 
paving the way towards sub-Nyquist Doppler focusing. 

Using Q, and denoting X v (uS) as the CTFT of x p {t+pr), 
we have 

L-x 

X p (w) = H{uj) ^2 aie- jUTl e- MpT . (9) 

1=0 

Taking the CTFT of (f>(t; v) as a function of t: 

p-i 

v) = CTFT($(<; v)) = ^ X v {u)é vvr 

p=0 

L-l P-l 

= H(u) Y meT 3UTl Y e j ^- v,)pr . 

1=0 p=0 



(10) 
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Ali 3L of the problem's parameters appear in (jTOj». Its structure 
is that of a delay estimation problem as we will see in 
Section Vi combined with the familiar sum of exponents term 
from |6i. 

We have seen that Doppler focusing reduces a delay- 
Doppler estimation problem to a delay-only problem for a 
specific Doppler frequency. In the next section we describe 
delay recovery from sub-Nyquist sampling rates using Xam- 



pling. In Section VI we revisit the Doppler focusing concept, 



combining it with these low rate Xamples to solve the delay- 
Doppler radar problem at sub-Nyquist rates. 

V. Sub-Nyquist Delay Recovery 
The problem of recovering the 2L amplitudes and delays in 



L-l 



(11) 



from sub-Nyquist samples has been previously studied in 0, 
0, 0, lfl2Ì . Since Doppler focusing yields such a problem, 
we now review how Xampling can be used to solve (jTTJ at a 
sub-Nyquist sampling rate. 

A. Xampling 

The concept of Xampling, introduced in m, lfl3ll and 
[14|, describes analog-to-digital conversion which acquires 
samples at sub-Nyquist rates while preserving the ability to 
perfectly reconstruct the signal. Xampling can be interpreted 
as "compressed sampling", in the sense that we are performing 
data compression inherently in the sampling stage. To do this, 
we do not simply reduce sampling rate, since this is bound 
to cause loss of information. Instead, we perform an analog 
prefiltering operation on our signal and only then sample it, in 
order to extract the required information for recovery. We now 
show how the signal's Fourier series coeffìcients are related 
to the problem's unknown parameters Q, 0, @, 1121 . We 
then describe how we acquire these Fourier coeffìcients via 
Xampling. 

Since (j>{t) is confìned to the interval t £ [0, r], it can be 
expressed by its Fourier series 



c[k]e 



j2-rrkt/T 



te [0,r], 



(12) 



where 



c[k] = - / 4>{t)e-^ kt/T dt 



l^ 1 



oti 



1=0 



: H(2irk/T 



h(t~ Tl )e-^ kt ^dt 



L-l 

aie J 11 . 

1=0 



(13) 



From ( fT3] l we see that the unknown parameters {a;, t;}^ 1 
are embodied in the Fourier coeffìcients c[k] in the form of 
a complex sinusoid problem. For these problems, if there is 
no noise, 2L samples are enough to recover the unknown a's 
and r's (9), i.e. |k| > 2L. We can solve this problem by 



using spectral analysis methods such as the annihilating filter 
ifBll or matrix pencil ED- One could also use MUSIC H3 
or ESPRIT fl8l . These recovery algori thms generally require 
sampling a consecutive subset of coeffìcients. The lower bound 
on |k| can be achieved only when the noise is negligible 
and computational complexity is not of concern. When there 
is substantial noise in the problem, having more than 2L 
coeffìcients will allow the recovery to be more robust. 

Our signals exist in the time domain, and therefore we 
do not have direct access to c[k]. We can use the Direct 
Multichannel Sampling scheme described in Q in order to 
obtain the Fourier series coeffìcients. Fig. [3] demonstrates 
a Xampling scheme used to directly extract the required 
Fourier coeffìcients from the signal. In an actual radar 
system was built using a similar yet more practical technique, 
where a set of mixers, band-pass filters and low rate ADCs 
sampled different spectral bands of the signal and the matching 
Fourier coeffìcients were created digitally. The radar model 
there included delay only without Doppler. An alternative 
Xampling method uses the Sum of Sincs filter as described 
in fl2l . Ali of these methods can be used to obtain arbitrary 
Fourier series coeffìcients. The number of Fourier coeffìcients 
extracted per pulse is a design parameter which controls the 
tradeoff between sampling rate and recovery performance. 
In our numerical experiments, presented in Section VII we 



demonstrate a Xampling rate one tenth of the Nyquist rate. 



-9— liC 



(:)dt 



(■)* 



Fig. 3: Multichannel direct sampling of the Fourier series 
coeffìcients, from 0. 



B. Compressed Sensing Recovery 

In the previous subsection we showed that 2L samples are 
enough to solve ( [11) , when there is no noise. We now describe 
a CS-based recovery method, operating on the Xamples c[k], 
which is more robust to noise. 

Assume the time delays are aligned to a grid 



ti = qiA T , < qi < N T 



(14) 



where we choose A r so that N T = r/A r is an integer. 
Choose a set of indices k = {fco, fci K i_i}, and defìne the 
corresponding vector of Fourier coeffìcients 



c = [c[*o]...c[/5 W _i]] T eCW. 



(15) 
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We can then write ( fT3j ) in vector forni as 
c = -HVx 



(16) 



where H is a |«| x |k| diagonal matrix with diagonal elements 
H(2irki/T) and V is a \k\ x N t matrix with elements ~V mq = 
e -j2irk m q/N T ^ j e j t j s com p 0se( j f | K | rows G f the N T x N T 

DFT matrix. The vector x e is L-sparse, where each 

index q contains the amplitude of a target with delay qA r if 
it exists, or zero otherwise. Defining the CS dictionary A = 
ÌHV e C^ xN - we obtain the CS equation 



c = Ax. 



(17) 



Estimating targets' delays can be carried out by solving ( fT7j ) 
and finding x's support - any nonzero index q denotes a target 
with delay qA T . 

For any set of sampled Fourier coefficients, a variety of CS 
techniques can be employed for recovery |10|, for instance 
orthogonal matching pursuit (OMP) 1191 . iterative hard thresh- 
olding (IHT) [20], or LI minimization dl2Tl and references 
within). Also, choosing the coefficients at random produces 
favorable conditions for CS, aiding recovery in the presence 
of noise. If the indices in k are selected uniformly at random, 
it can be shown that if |k| > CL(logN T ) 4 , for some positive 
Constant C, then A obeys the desired Restricted Isometry 
Property (RIP) with large probabili ty \22\. By satisfying the 
condition for RIP we are able to recover x, using a CS 
recovery algorithm. 



VI. 



Delay-Doppler Recovery Usino Doppler 
Focusing 



In Section IV we introduced the concept of Doppler fo- 



cusing, and in Section [V] we reviewed how to Xample and 
recover the 2L unknowns of the delay-only problem ( fTTj i. We 
now return to our originai delay-Doppler problem |2|. 

We begin by describing how the target scene can be re- 
covered perfectly when there is no noise, at a sample rate 
proportional to the signal's rate of innovation. This recovery 
method is not robust to noise, motivating the use of Doppler 
focusing, described after we show how Xampling can be 
performed on the multi-pulse signal Q. Finally we discuss 
some practical considerations and compare our method to 
previous CS approaches. 

A. Noìseless Recovery 

We are interested in analyzing the noiseless scenario in 
order to find the minimal number of samples K m i n needed 
for perfect recovery of our problem' s 3L unknown parameters. 
Since K m i n cannot be smaller than the number of DOF, 
K m in > 3L. We now describe a reconstruction method 
capable of perfectly estimating ali 3L unknowns using 4L 
samples, slightly above the minimal rate but stili significantly 
below the Nyquist rate for reasonable £?/, and L. 

Defining /3f = aie~i" lVT we obtain 



L-l 



c p [k] = -Hpnk/T^fté 



-j2irkT l /T 



(18) 



Choosing a set of indices k for which H(2nk/T) ^ 0, we can 
preprocess the coefficients c p [k) by the inverse i7^ 1 (27rfc/r) 
to obtain a standard spectral analysis problem, as we have seen 
in Section [V] For these problems, 2L samples are enough to 
perfectly recover the unknown /3's and r's, if there is no noise. 

Suppose we have used 2L samples to solve ( fT8j ) for some 
< p < P — 2. We now posses the set of delays {t;}^, 1 , thus 
recovering the first L DOF. We also posses the set {j3f }f=Q, 
to be used shortly. To recover the remaining 2L DOF, we solve 
( fT8| ) for p' =p+ 1, obtaining {/Sf }^ 1 . Since 



aie 



-jvi (p+l)r 



JUIT 



we can recover the set of Doppler frequencies by 



Vi 




(19) 



(20) 



Using ( |20] >, recovery of the set of amplitudes is immediate 
from ai = ffie 1 "^. 

To conclude, we have shown a reconstruction method for 
the noiseless scenario, which perfectly recovers the problem's 
3L unknown parameters using only AL samples. In order to 
provide robustness to noise, we now combine Xampling with 
Doppler focusing. 



B. Xampling 

Similarly to the Xampling technique of Section [V] which 
obtained c[k], we can extend this technique to each of the 
pulses x p (t) of the multi-pulse signal |2]) to obtain c p [k]. Since 
x p {t) is confined to the interval t E [pr, (p + l)r], we can 
replace t — > t + pr and ai — > aie~ì" lpT in ( fT3| l to obtain 



L-l 



Cp[k] = -H(2nk/T) V aie- jVlpT e- j2wkTl/T , 



2=0 



(21) 



where we used the fact that since both k,p E Z we have 
e -j2Trkp = ^ From pTj ) we see that ali 3L unknown param- 
eters {ai, Ti, vi}^Sq are embodied in the Fourier coefficients 
c p [k) in the form of a complex sinusoid problem. 

C. Applying Doppler Focusing and CS Recovery 

Having acquired c p [k] using Xampling, we now perform the 
Doppler focusing operation for a specific frequency v 

p-i 

¥„[*] - c vW 3VVT 



p=0 



L-l 



p-1 



= -H{2nk/T) J2 aie^ kTl/T é {v ~ Vl)vT . (22) 

1=0 p=0 

From ( fT0] > we see that ^ y [k] = ^(u>; v)\ LU =2Trk/T- 

Following the same arguments as in (|7), for any target l 
satisfying \v — vi\ < 2it/Pt we have 



p-i 

p=0 



(23) 
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Therefore, Doppler focusing can be performed on the low rate 
sub-Nyquist samples: 

V v [k] S -H(2irk/T) V ai e- j2vkT '/ T . (24) 

i:|«/-^i|<27r/Pr 

Equation ( p4| i is identical in forni to ( pj) except it is scaled 
by P, increasing SNR for improved performance with noise. 
Furthermore, we reduced the number of active delays. For 
each v we now have a delay estimation problem, which can 
be written in vector forni using the sanie notations of Section [V| 
as 

*„ = -HVx„ (25) 

T 

where 

= ...^[fc w _ifeCW. ( 26) 

This is exactly the CS problem we have already shown how 
to solve, with one important difference. In the delay-only 
problem of Section [V] the model order L was known. Here, 
since there are L targets but we have no prior knowledge of 
their distribution in the delay-Doppler grid, for each v we 
must either estimate the model order < L v < L, or take 
a worst case approach and assume L v = L. The problem of 
estimating the number of sinusoids in a noisy sequence has 
been studied extensively ||23Ì . Il24l . Il25ll . Solving ( |25] l with an 
accurate model order can decrease computation time (although 
estimating model order is also time consuming) and possibly 
reduce detection of spurious targets. The former option is 
preferable in higher SNR scenarios when the L v estimates 
are good, while the latter can be a fall back approach in noisy 
scenarios when the L u estimates contain signifìcant error. In 
our simulations, since we do not want model order errors to 
influence recovery performance, we employ the worst case 
approach. 

The Doppler focusing technique is a continuous operation 
on the variable v, and can be performed for any Doppler 
frequency up to the PRF Since the focus zone for each v 
is of width 2-k/Pt, we can fìnd various finite sets of i/'s 
spanning [0,2tt/t]. For any such set, defìne its size as N v . 
For each v in the set, we solve ( |2"5| > assuming x„'s support 
is of size L. This problem can be solved using an abundance 
of CS algorithms as described in Section [V] After solving N v 
separate CS problems with dictionary of size \n\ x N T , we 
hold at most LN V estimated amplitudes. Since the absolute 
value of amplitudes recovered in the support is indicative of 
true target existence as opposed to noise, we sort the LN„ 
estimated amplitudes according to their absolute value, and 
take the L strongest ones as true target locati ons. Algorithm [T] 
summarizes the Doppler focusing algorithm. 

D. Practical Considerations 

If the set of probed Doppler frequencies lies on a uniform 
grid, i.e. v n = 2im/N v T,n = 0,1,..., N v - 1, then 9 v [k] 
can be created effìciently using a length N v DFT or FFT of a 



Algorithm 1 Doppler Focusing 

Input: Xamples {c p [k]}^ <p 

Output: Estimated target parameters {ài, 77, z^}^ 1 

Initialization: Discretize unambiguous Doppler interval 
[0, 27t/t] into N v points {^nj^Q 1 , where the union 
of focus zones of ali points covers the interval. 

for n = to N v - 1 do 

Create \I>„ n from {cp[fc]}o<p<p using Doppler focusing 
(f22| 

Solve ^ Vn = ^HVx„ for x„ to obtain {ài,fp}f~J 
end for 

Sort {|à"|}g!™ < f^ V ", denoting indices of L largest values 

as {li, rii}f_Q . 
for i = to L — 1 do 

ài <- à™. 1 ; n <- ff* ; h v n% 
end for 



series of length P: 

p-i 

p-i 

= ^ c p [k]e^' N ^ = DFT Nv {c p [k}}. (27) 

p=0 

When focusing on some Doppler frequency v, targets with 
Doppler frequencies v\ satisfying \v\ — v\ > 2tt/Pt are 
considered undesirable and we seek to minimize their effect. 
These targets can be viewed as "out-of-focus", since they are 
not matched to v and their responses from different pulses 
do not combine coherently; they will combine in phase for 
different v's satisfying \vi — v\ < 2tt/Pt. We can add to p2| 
a user defìned weighting function w[p],p = 0,1,..., P — 1 
(e.g. Hann, Blackman, etc.) which is designed to mitigate the 
impact of these out-of-focus targets: 

p-i 

*u[k] = cp[k]^w\p] 

L-l P-l 

= -H(2irk/T) ^2 aie- j2wkTl / T ^ e 3{ "~ Vl)pT w[p\. (28) 
T 1=0 p=0 

The drawback of weighting is that it increases the frequency's 
focus zone, potentially including more targets in each delay 
estimation problem. In Fig. |4] we see an example of how 
weighting functions can reduce the effect of out-of-focus 
targets compared with no weighting (Constant For a 

comprehensive review of weighting function design consider- 
ations see [26|. 

Slow time undersampling, i.e. sampling only a subset of 
the P transmitted pulses, is a naturai extension of the Doppler 
focusing method. Instead of Xampling ali P pulses to obtain 
c p [k] for p — 0,1, ...,P — 1, and then combining them ali 
as in ( |22] i, we can work with a subset of the pulses T C 
{0, 1, P — 1}. Undersampling using T is advantageous in 
the same ways as fast time undersampling: reducing hardware, 
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- Constant w[p] 

- Hann 

- Blackman 




-0.2 -0.1 0.1 0.2 

Normalized Doppler frequency v 

Fig. 4: DFT of weighting functions w[p] compared with no 
weighting (Constant w[p\) for P = 100 pulses. Attenuation 
of targets with Doppler frequencies far enough from nominai 
frequency (v = in this case) increases significantly with 
proper weighting. Also, focus zone changes for different 
u>[p]'s. 



processing and Storage requirements. The disadvantage of slow 
time undersampling is also like that of fast time - a decrease 
in SNR and possible decline in estimation performance. It 
adds another degree of freedom to the sampling and recovery 
architecture, allowing to split the undersampling ratio between 
fast and slow time samples and balancing the tradeoff between 
them. For example, a system with Storage constraints, may use 
F to dilute the number of processed pulses while maintaining 
a higher fast time s ampie rate. 

E. Prevìous CS Approaches 

A possible approach to the problem at hand could be a two- 
stage CS recovery technique, first estimating target delays and 
afterwards for each delay estimating Doppler frequency (or 
vice versa). These two-stage methods tend to be suboptimal 
since the problem is not separable, and a mistake in the first 
estimation stage propagates to the second stage where it cannot 
be undone. We compare our method to this type of recovery 
in Section IVIII 

To overcome this inefficiency, several works 0, O employ 
a single stage CS recovery technique. Instead of estimating de- 
lays and Doppler frequencies separately and sequentially, they 
estimate the most likely (17,2/;) pairs. The drawback of this 
technique is that it requires using a dictionary with dimensions 
proportional to N T N V F S , where F s is the problem's sample 
rate. Since grid sizes can easily reach IO 3 , and sample rates are 
on the order of MHz, the dictionary grows rapidly rendering 
these methods infeasible for even moderate problem size. 

VII. SlMULATION RESULTS 

In the previous sections we described the Doppler fo- 
cusing approach to delay-Doppler recovery. Before showing 



numerical examples comparing our method to other recovery 
techniques, we first discuss how the user defined performance 
metric infiuences grid size and Fourier coefficient selection. 

A. Performance Metric 

Our problem lies in a continuous, analog world. Since we 
choose to solve it using CS, which is an approach developed 
for discrete problems, we must discretize the delay grid, to 
steps of size A T . As real world targets delays and Doppler 
frequencies do not He on any predefined grid, but our CS 
recovery assumes they do, it seems we should take A T — > 
in order to minimize discretization errors. Computational 
requirements aside, there is one significant drawback to such a 
decrease in grid step - columns in the CS dictionary A become 
increasingly similar, making A increasingly coherent, where 
coherence is defined as the largest absolute inner product 
between any two columns aj, aj of A: 

Ka,,a,)| 



M(A) 



max t 

|ai||2||aj|| 2 



(29) 



A basic premise of CS ties low coherence to successful 
recovery [27 1. Therefore, A is usually designed to have small 
coherence. This contradicts taking the step size to be very 
small. 

Here, we move away from this basic assumption, and 
argue that depending on the chosen performance metric, high 
coherence can actually help recovery instead of harm it. For 
example, assume we are interested in delay recovery but are 
tolerant of some small error r max . In radar applications, a 
common performance metric is the "hit-or-miss" criterion on 
the estimated delays {t;}^, 1 : 



if \n - n\ < r max 

otherwise 



L-l 

E 

1=0 



et. 



(30) 



Translating r max to a condition on support recovery, ( |30| ) 
tolerates an error of no more than K — \r max / A T \ E N 
places in the recovered indices from (17) , where we assume 
K > 1. Instead of designing A so that each column is as non- 
correlative with the other columns as possible, we can try to 
design the dictionary so that each column is somewhat correl- 
ative with its K nearest neighbours, and only afterwards does 
the correlation drop. This will improve recovery performance 
in noisy scenarios since in cases where the correct column is 
not recovered, any one of its similar neighbouring columns 
stili has a chance to overcome the noise and produce a "hit". 
Graphs in the next subsection will show how such a coherent 
A actually improves recovery performance compared with a 
less coherent dictionary in very noisy scenarios. 

If T max is chosen so that K = 0, the same line of thought 
entails requiring each dictionary column to be similar to 
neighbours, i.e. for A to have minimal coherence. Since most 
CS works deal with exact recovery, this explains why they 
strive for a minimal p(A). 

We can control the level of A's coherence by choosing 
different sets of Fourier coefficient k in ( fl5] l. Fig. [5] shows 
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an example of the column correlation pattern for two sets of 
Fourier coefficients: a consecutive set and a random set, where 
ali coefficients were chosen in [—B^/2, B^/2}. We define the 
column correlation function for some column ai as 



l( a i, a j)l 

I ai 1 1 2 1 1 aj 1 1 2 ' 



(31) 



The consecutive set is better suited for performance criteria 
which allow some error in support recovery, while the random 
set will achieve better performance when exact recovery is 
required. 




-Consecutive coefficients 
- Random coefficiente 




yvwyj 



Fig. 5: Column correlation pattern no[j} of the CS dictionary 
from ( fT7| for two sets k. The consecutive set achieves coher- 
ence of 0.9 due to the many correlative columns in the center, 
while the random set has coherence of 0.3. 



method as described in [7| (where we use a CS algori thm 
instead of ESPRIT) using the following criteria: 

1) Hit-Rate - we define a "hit" as a delay-Doppler estimate 
which is circumscribed by an ellipse around the trae 
target position in the time-frequency piane. We used an 
ellipse with axes equivalent to ±3 times the time and 
frequency Nyquist bins. 

2) Recovery RMS error - for estimates classified as hits, 
we measure the root mean square error in both time and 
frequency. 

As noted in the previous section, a single stage CS recovery 
method using Nyquist bins spacing consumes a prohibitive 
amount of memory and was not able to run on any computer at 
hand, since the CS dictionary required storing 4- IO 9 elements 
(occupying 32GB of memory using standard IEEE doublé 
precision): ^j/pr i/b, = ^ >T ^h = 2 ■ IO 5 columns and 
PrBh/ '10 = 2 ■ IO 4 measurements per column. 




B. Numericeli Results 

We now present some numerical experiments illustrating the 
recovery performance of a sparse target scene. We corrupt the 
received signal x(t) with an additive white Gaussian noise 
n(t) with power spectral density S n (f) = No/2, bandlimited 
to x(t)'s bandwidth B^. We define the signal to noise power 
ratio for target l as 



SNR, 



^f^\ ai h(t)\ 2 dt 



(32) 



N B h 

where T. p is the pulse time. The scenario parameters used 
were number of targets L=5, number of pulses P=100, PRI 
r=10/isec, and B^=200MHz. The classic time and frequency 
resolutions ("Nyquist bins"), defined as l/B^ and 1/Pr, are 
5rtsec and 1 KHz accordingly. In order to demonstrate a 1:10 
sampling rate reduction, our sub-Nyquist Xampling scheme 
generated 200 Fourier coefficients per pulse, as opposed to 
the 2000 Nyquist rate samples. We tested Doppler focusing 
with two types of Fourier coefficient sets k, a consecutive set 
and a random set. We compared Doppler focusing recovery 
performance with classic processing and a two-stage recovery 



Fig. 6: Hit Rate for classic processing, two-stage CS recov- 
ery and Doppler focusing as function of SNR. Sub-Nyquist 
sampling rate was one tenth the Nyquist rate. 

For CS-based techniques, the delay grid step A T was chosen 
as half a Nyquist bin. For Doppler focusing, the Doppler 
frequency region was discretized with uniform steps of half a 
Nyquist bin. To provide a fair comparison, classic processing 
performed the slow time FFT after padding with P zeros. 

Fig. [6] and Fig. [7] demonstrate the hit-rate and RMS error 
performance of the different recovery methods for various 
SNR values. It is evident that Doppler focusing is superior 
to the other sub-Nyquist recovery techniques. Between the 
two Doppler focusing approaches, consecutive coefficients 
are better suited for lower SNR, while choosing coefficients 
randomly improves performance as SNR increases. Since both 
sets k produce CS dictionaries with column correlation func- 
tions which are not matched to the "hit-or-miss" performance 
criteria used, we have no reason to assume one should be better 
than the other. Also, random coefficients, when producing a 
hit, have very small delay errors (even compared with Nyquist 
rate classic processing) due to low CS dictionary coherence. 



IO 



Time RMS errar 




; rjrocesein B (Nvquist) 
Classic processing(Sub-NyquistJ 
Two stage CS(Sub-Nyq) 
Doppler Focusing - consecutive coeff.(Sub-Nyq) 
Doppler Focusing - rsndom coeff (Sub-Nyq) 




Fig. 7: RMS error of time and frequency estimates for classic 
processing, two-stage CS recovery and Doppler focusing as 
function of SNR. Sub-Nyquist sampling rate was one tenth 
the Nyquist rate. 

Fig. [8] shows the sparse target scene on a time-frequency 
map for a -28dB SNR scenario, where each target is displayed 
along with its hit rate ellipse, together with the various sub- 
Nyquist recovery methods' estimates and hit rates. As noted 
only Doppler focusing is able to distinguish 



in Section IV 



between the two targets having almost identical delays (around 
4.2 /xsec) but different Doppler frequencies. 

L=5 P=100 Avg Sample SNR=-2B D[dB] 
Hit Rate=D 6D(MF-SubNyq) 20[Two stage CS-SubNyq) 0.80 (Doppler focusing-SubNyq) 



Real 

Classic processing (Sub-Nyquist) 
Two-stage CS (Sub-Nyquist) 
Doppler Focusiog (Sub-Nyquist) 



5 6 
time [micro-secj 



Fig. 8: Real target positions along with various estimates. 
Doppler focusing achieves highest hit rate among sub-Nyquist 
methods. Only Doppler focusing detects two targets around 
4.2/isec. 



Vili. CONCLUSION 

We demonstrated a radar sampling and recovery method 
called Doppler focusing, which employs the techniques of 



Xampling and CS, and is independent of the radar signal's 
bandwidth. Doppler focusing allows for low rate sampling 
and digitai processing and undersampling in both fast and 
slow time. It also leads to CS recovery with dictionary size 
scaling with delay grid size only, and imposes no constraints 
on transmitted signal or modulation. We compared our method 
to other sub-Nyquist recovery techniques and have seen its 
clear advantage in simulations. When sampling at one tenth 
the Nyquist rate, and for SNR above -25dB, Doppler focusing 
achieves results almost identical to classic recovery working 
at the Nyquist rate. 

Appendix 

In the following section we formally define assumptions Al- 
A3. For this we define a few auxiliary parameters that were 
previously unneeded. Denote the radar' s carrier frequency as 
/ c , the pulse time support as T p , and the speed of light as 
c. We use the time-distance equivalence ri — cti/2 and the 
non-relativistic Doppler radiai frequency-velocity equivalence 
ri = uic/ATrf c . 

Al. "Far targets" - target-radar distance is large compared to 
the distance change during observation interval 



riPr < r t 
vicPt/AkJc < cti/2 

vi < 2nf c Ti/PT. 



(33) 



A2. "Slow targets" - we can assume Constant ti if the base- 
band Doppler frequency is smaller than the frequency 
resolution 



2riB h /c « 1/Pt 
viB h /2irf c « 1/Pt 

vi « 2-Kf c /PrB h . 



(34) 



A3. 



This inequality is termed the "narrowband assumption" 
in radar nomenclature, since it gives an upper limit on 
the radar signal's bandwidth. 

We can assume the Doppler induced phase change over 
a single pulse time is small if v{F v <C 1. 
"Small acceleration" - velocity change induced by ac- 
celeration is smaller than velocity resolution 



nPr « c/2f c PT 

h « c/2f c (Pr) 2 . (35) 
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